STIRAP transport of Bose-Einstein condensate in triple-well trap 
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The irreversible transport of multi-component Bose-Einstein condensate (BEC) is investigated 
within the Stimulated Adiabatic Raman Passage (STIRAP) scheme. A general formalism for a 
single BEC in M-well trap is derived and analogy between multi-photon and tunneling processes 
is demonstrated. STIRAP transport of BEC in a cyclic triple-well trap is explored for various 
values of detuning and interaction between BEC atoms. It is shown that STIRAP provides a 
complete population transfer at zero detuning and interaction and persists at their modest values. 
The detuning is found not to be obligatory. The possibility of non-adiabatic transport with intuitive 
order of couplings is demonstrated. Evolution of the condensate phases and generation of dynamical 
and geometric phases are inspected. It is shown that STIRAP allows to generate the unconventional 
geometrical phase which is now of a keen interest in quantum computing. 

PACS numbers: 03.75.Mn, 03.75.Lm, 05.60.Gg 



I. INTRODUCTION 

In recent decades, investigation of Bose-Einstein con- 
densate (BEC) has become one of the main streams 
in modern physics (see for reviews and monographs 
P, 0, H, 0) H, H 0])- I n particular, a large attention is 
paid to dynamics of multi-component BEC and related 
coherent phenomena, e.g the Josephson effect in weakly 
bound BECs . 

Two kinds of multi-component BEC are usually con- 
sidered. The first one is confined in a single trap and 
contains atoms in a few hyperfine levels (multi-level sys- 
tem or MLS). Here every component is formed by atoms 
in a given level. The components can be coupled by the 
laser light with the carrier frequency equal or close to 
the difference of the Bohr frequencies of the hyperfine 
states. One can control the coupling by varying parame- 
ters of the laser irradiation and so get different regimes of 
the transfer of atoms between the components: Joseph- 
son oscillations (JO), macroscopic quantum self-trapping 
(MQST), etc (see e.g. [lOj for discussion). 

In the second kind of the multi-component BEC, the 
atoms are in the same hyperfine state but the trap is sep- 
arated by laser-produced barriers into a series of weakly 
bound wells (multi-well system or MWS) [1J]. BEC 
atoms can tunnel through the barriers and exhibit the 
similar effects as the MLS. In this case BEC components 
are represented by populations of the wells. Both MLS 
and MWS are obviously similar. Indeed, the coupling 
Rabi frequencies in MLS are counterparts of the barrier 
transition matrix elements in MWS. And detuning be- 
tween Bohr and carrier frequencies in MLS is similar to 
detunin g of the well depths in MWS. BEC in optical lat- 
tice (25l |26| can be also treated as MWS though, unlike 
a few- well traps [24J, the well depths and separations in 
optical lattices cannot be usually monitored individually 
for every cell. 

Most of the studies consider BEC with two components 
[8||-[ll| and much less with three components fl6|-[23j|. 



The later case is much more complicated. At the same 
time, it promises new dynamical regimes (l7l . [Lsj and 
allows to consider not only linear (couplings 1-2, 2-3) 
but also cyclic (couplings 1-2, 2-3, 3-1) well chains. 

In the present paper, we investigate BEC dynamics in 
the triple-well trap, i.e. MWS with the number of wells 
M=3. Unlike most of the previous studies, we will ex- 
plore not oscillating fluxes of BEC but its complete and 
irreversible transport between the initial and target wells. 
For this aim, the coupling between BEC fractions ^com- 
ponents) will be monitored in time (unlike the constant 
coupling for the Josephson-like oscillations). Once being 
realized, BEC transport could open interesting perspec- 
tives in many areas, e.g. in exploration of coherent topo- 
logical modes [13, HI] and diverse geometric phases (29j - 
|32| . The later is especially important since geometric 
phases are considered as promising information carriers 
in quantum computing [33l |34| . l35j . 

Due to similarity between MLS and MWS, one may 
try to apply for BEC transport numerous develop- 
ments of modern quantum optics, in particular, adiabatic 
two-photon population transfer methods 36]. Between 
them the stimulated Raman adiabatic passage (STIRAP) 
[361 . H3] seems to be the most suitable for our aims since 
it allows, at least in principle, the complete irreversible 
population transfer. The method was first developed for 
atoms and simple molecules [36l . W\ and then probed 
in metal clusters l38| -[4fj| and variety of other systems, 
see references in |22|. Quite recently STIRAP was ap- 
plied to the tr ansp ort of individual atoms [12] and BEC 
[U, El [U, [H, HI in the triple- well trap. 

The applicability of STIRAP to BEC transport needs 
some care since interaction between BEC atoms results 
in a time-dependent nonlincarity of the problem, which 
can destroy the adiabatic transfer 0, Hi], H2, HH ■ This 
nonlincarity plays the same detrimental role as the dy- 
namical Stark shift in electronic MLS systems, where it 
disturbs the two-photon resonance condition and thus 
breaks one of the basic STIRAP requirements (see dis- 



cussion in Sec. |TT|. As was shown in (l9| . the undes 
nonlinear impact can be circumvented by using a ( 
ing larger than the atomic interaction. The subse 
studies [22, 23] confirmed that the detuning is usefu 
aim a minimal (say < 1%) occupation of the interm 
well during STIRAP process. The less the (tempi 
occupation, the better adiabaticity and robustness 
process. At the same time, one should recognize th 
cupation of the intermediate well cannot be fully av 
Moreover, that occupation is temporal and in any i 
further transferred to the final well. So it does not 
the final fidelity of the BEC transport. 

In this study we will show that the robust and cor 
transport of the interacting BEC can take place e' 
zero detuning, regardless of the temporary weak pi 
tion of the intermediate state. Of course such trai 
is more likely quasiadiabatic but we are interested 
transport completeness rather than in its purely 
batic character. Moreover, we will show that the 
plctc transfer can be done even at intuitive seque 
the pump and Stokes couplings (unlike their cour 
tuitive order in STIRAP), i.e. in strictly non-adi 
case. 

As compared with the previous studies [Til HT| Q 
we will consider more general triple-well trap whi< 
also 3-1 cou plin g and thus represents the circular c 
uration [ltjl |21|. Such configuration allows to run 
through the circle as many rounds as we want ar 
it to any of three wells. Besides the population 
condensate phases will be explored. Moreover, v, 
present some first examples of the dynamical ana geo- 
metric phases generated in STIRAP. The later is possi- 
ble because the circular well configuration and adiabatic 
STIRAP transfer allow to build the adiabatic cyclic evo- 
lution. It worth noting that condensate phases and their 
dynamical and geometric constituents were not yet ex- 
plored in STIRAP (for exception of a brief phase analysis 
in [U). 

The paper is outlined as follows. In Sec. II we sketch 
STIRAP. In Sec. Ill a general mean-field formalism 
for dynamics of multi-component BEC is presented and 
specified for MWS with M=3. In Sec. IV the calculation 
scheme is given and similarity between our scenario and 
conventional STIRAP is discussed. In Sec. V results of 
the calculations are discussed. The conclusions are done 
in Sec. VI. 



II. STIRAP 

STIRAP [H, H?! is the adiabatic two-step process pro- 
viding the complete population transfer from the initial 
level |1) to the target level |3) via the intermediate level 
1 2). Its scheme for MLS is illustrated in Fig. [TJ As 
is seen, the transfer is driven by the pump and Stokes 
laser pulses with Rabi frequencies flp(t) and fls(t)- The 
pump laser couples the initial and intermediate states 
while the Stokes laser stimulates the emission from the 



a) l 2 > — >A b) 





FIG. 1: a) STIRAP scheme for a three-level A-system. The 
pump and Stokes pulses with Rabi frequencies Qp(t) and 
Sls(t) provide the couplings |1) — 12) and |2) — 13), respectively. 
A is the detuning from the intermediate level, b) Counterin- 
tuitive sequence of pump and Stokes pulses overlapping for a 
time t. 



intermediate state to the target one. In addition to the 
A-configuration given in the plot a), STIRAP can be also 
realized in the ladder and ^-configurations [36| . 

STIRAP has three principle requirements: 
1) two-resonance condition for the laser carrier and Bohr 
frequencies 



LUp — U>S = — U>l 



(1) 



allowing a detuning A = (lu2 — a>i) — ujp = (uj2 — ^3) — lus 
from LU2] 

2) overlap of the pulses (during the time r) and their 
counterintuitive order (the Stokes pulse proceeds the 
pump one); 

3) adiabatic evolution ensured by the condition 



VLt > 10 



(2) 



where f2 = y/flp + Q s . 

Due to interaction with the laser irradiation, the bare 
states |1), 1 2), |3) are transformed to the dressed states 

\a + ) = sin 9 sin (f>\ 1) + cos 2) + cos sin cf>\ 3), 
\a°) = cos0|l) -sin0|3), (3) 
|a~) = sin#cos</)|l) — sin0|2) + cos#cos</>|3) 



with the spectra 



= A± y/A 2 + n 2 , 







(4) 



and mixing angle determined through 

sine = n P /ri, cosd = ri s /n. (5) 

The mixing angle <fr, a known function of the Rabi fre- 
quencies and detuning [41'], is of no relevance in the 
present discussion. For the sake of simplicity, we omitted 
in this section time dependence of Rabi frequencies and 
other values. 

STIRAP has been developed for population of non- 
dipole states (with the spin J / 1) which cannot be 
excited in the photoabsorption but can be reached by 
two dipole transitions via an intermediate state. The 
main aim was to avoid, as much as possible, the pop- 
ulation of the intermediate state and thus prevent the 
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leaking via decay of this state. As it is seen from Eq. 
©, this can be achieved by keeping the system dur- 
ing all the process in the dressed state |o°) which has 
no contribution from the intermediate state |2). For 
this aim we need an adiabatic evolution following with 
the counterintuitive sequence of the pump and Stokes 
pulses as is shown in Fig. lb. In this case, we have 
^s(t) 7^ 0, ttp(t) = 0,sin# = 0,cos# = 1 at early time 
and fl s (t) = 0,fl P (t) ^ O,sin0 = l,cos0 = at late 
time. Then, unlike \a + ) and \a~), the state \a°) is re- 
duced to |1) at the beginning and to |3) at the end of 
the evolution. The main point is to evolve the system 
adiabatically, keeping it all the time in the state |ao). As 
is seen from ([2]) , for this aim we need either a strong cou- 
pling O (i.e. high laser intensities) or a long overlapping 
time r which is just duration of the adiabatic evolution. 
The STIRAP Hamiltonian reads 

n ( o n P 0\ 

H(t) = - [ftp 2A Q s . (6) 

1 \ o n s o J 

Note that eqs. ©-(jH) are obtained in the rotating wave 
approximation (RWA). Using the fact that u>p t s, ^1,2,3 3> 
Qp,s, A, the RWA allows to omit the high laser and Bohr 
frequencies and keep only the low frequencies of interest. 

Eqs. ©-© also neglect the dynamical Stark shifts 
6i(t) pertinent to MLS, which are supposed here to be 
weak, Si(t) <C flp t s, A. Otherwise the shifts enter the di- 
agonal terms in ^ and complicate Eqs. AS])-© ,36]. The 
dynamical Stark shifts are detrimental for STIRAP since 
they destroy the two-resonance condition ((T|) . The larger 
Rabi frequencies flp y s (and hence the laser intensities), 
the stronger the Stark shifts. So, £lp ; s must be large 
enough to keep the adiabatic condition and, at the 
same time, small enough not to cause too strong Stark 
shifts. This problem is obviously absent in MWS. As 
was mentioned above, the time-dependent nonlinearity 
in BEC is detrimental for STIRAP 0. In fact it plays 
the similar destructive role as the Stark shift. 



III. MEAN-FIELD DESCRIPTION OF BEC 
DYNAMICS IN MWS 

We start from the non-linear Schrodinger, or Gross- 
Pitaevskii equation [43] 

ih*(r,t) = {-^V 2 +V ext (r,t)+g \*(r,t)\ 2 Mr,t) (7) 
2m 

where the dot means time derivative, ^f(r,i) is the clas- 
sical order parameter of the system, V ex t (r, t) is the 
external trap potential involving both (generally time- 
dependent) confinement and coupling, go = Ana/m is 
the parameter of interaction between BEC atoms, a is 
the scattering length and m is the atomic mass. 

In what follows, we will consider the MWS where BEC 
is distributed between M wells separated by barriers. 



Then BEC components are reduced to the condensate 
fractions (=populations) in the wells. For BEC with 
weakly bound M fractions, the order parameter can be 
expanded as 

M 

#(f,t) = \/iV^ fe (t)$ fe (f) (8) 

k=l 

where $fc(r) is the static ground state solution of ([7]) for 
the isolated (without coupling) k-th well [3] and 

Mt) = v/A^iJe^W (9) 

is the amplitude related with the relative population 
Nk(t) and corresponding phase 4>k(t) of the k-th frac- 
tion. The total number of atoms N is conserved: 

jdWf,t)\yN = Y:tiN k (t) = i. _ 

Being mainly interested in evolution of populations 
Nk(i) and phases (j>k{t), we dispose by integration of the 
spatial distributions $fc(f) and finally get |9l.[To| 

M 

= [E k {t) + C/AI^I 2 ]^ - J2 Q kj{t)^ (10) 

where 

n kj (t) = -ij dr ■ v$i + nv**t(t)*i] (11) 

is the coupling between BEC fractions, 

Ek(t) = \l dr[^-\Vn\ 2 + <!>lV ext (t)<S> k } (12) 
a J 2m 

is the potential depth, and 

U=^Jdr\<P k \ 4 (13) 

labels the interaction between BEC atoms. The values 
f2fc,-(i), Ek(t), and U have dimension of the frequency. 

For simplicity we suppose that all the couplings have 
a common peak amplitude K. Then it is convenient to 
pick out this amplitude from the couplings 

fl kj (t) = KCl kj (t) (14) 

and scale (ITU1) by l/2if so as to get 

1 M 

ii> k = [E k {t) + A|^ fe | 2 ]Vfc - - J2^j(t)ipj • (15) 
Here 

E k {t) = E k {t)/2K, A - UN/2K (16) 

and the time is scaled as 2Kt — > t thus becoming dimen- 
sionless. Eq. (|15p is convenient since it is driven by one 
key parameter A responsible for the interplay between 



the coupling and interaction. As is shown below, th 
parameter is decisive for STIRAP transport of BEC. 

By substituting J9]) into (fl5|) and separating real an 
imaginary part, one gets equations describing evolutio 
of the system in terms of fractional populations N k (~i 
and phases (f> k (t): 



M 



N k = - ^ ttkjy/NjNk sm(<^- - 4> k ) 



(1< 



4>k 



1 M 




--ki\l ~rr~ COS( 
2 ^ 3 V N k V 



In MWS, the condensate is distributed between I 
space-shifted coupled wells, Vt r a P (r) — > Y^k=i Vk(r), wit 
the depths E k . The wells are separated by the barriei 
with penetrabilities fl k7 tj(t). We consider below a wea 
coupling between BEC fractions. Then only the neighbc 
fractions are coupled, Cl k j ^ for 1 < j = A; ± 1 < h 
and interaction between atoms of different fractions ca 
be neglected. 

Equations (fT7l) - (fT8| allow the classical analogy with 
the populations N k (t) and phases 4> k (t) treated as clas- 
sical canonical conjugates. It is easy to verify that these 
equations can be casted in the canonical form 



ggcj 

d<t> k 



(19) 
(20) 



dN k 

with the classical Hamiltonian 

M A/V 2 
H d = - Y J {EkN k + -^-) 

k 

+ 2 X) VkjVNkNj cos(& - ct> k ) . (21) 
kj 

One may further upgrade (fl~7|) - (|18p by means of canon- 
ical transformation of N k and cj) k to canonical pairs re- 
lated to the population imbalances and phase differences. 
This will allow to remove from (|17p - ([T8"l) the integral of 
motion N and decrease the total number of equations 
from 2M to 2(M-1). 

Let us consider the linear canonical transformation 



Z k 



kj 



(22) 



with Zk and 6 k being the population imbalances and 
phase differences, respectively. For the linear transfor- 
mation, the matrices T and R in (f2"2"| are related as 



R = T- 



(23) 




© © © 



FIG. 2: a) The linear triple-well system containing the in- 
teracting BEC. The wells are denoted by the numbered cir- 
cles. The tunneling couplings Q.vi(t) and Q.23{t) are marked 
by double arrows. Ei , E2 and E3 are well depths, b) The cir- 
cular configuration of triple-well system with the additional 
coupling 1-3. The solid lines connecting the circles mean the 
couplings. The dash curves represent the successive STIRAP 
transfers (1), (2) and (3), after which BEC returns back to 
the initial well |1). At every transfer only two couplings are 
active, as is shown in the plot a). 



For the MWS case with M=3 it is natural to chose 



zi = N 2 - 


Ni, 


Z2 


= N 3 - 


N 2l 


Z3 = 


N . 


Then one gets 
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(24) 



(25) 



Matrices for the inverse transformation 



M 



for M=3 read 

T -i = 



E R kj ld 3 



(26) 





-1 













■ 


-1 


; 






-1 


i 




2 






V 


1 





(27) 



New variables include the integral of motion z 3 = N 
and the total phase $3 = cj>i + 4>2 + <f>z. The equations for 
these values are straightforwardly separated from (fl7|) - 
(fT8|l and can be skipped. Then it is enough to solve the 
remaining four equations for z\, Z2, Q\ and 62 . 

Note that the formalism presented above is general and 
can be applied to both i) oscillatingJSEC fluxes in traps 
with constant parameters (like in [H"[l5| and [l7l|- [20l| ) 
and ii) irreversible BEC transport in traps with time- 
dependent parameters, e.g^ p, k j(t). This formalism is 

partly given elsewhere [1]-[23|]. However, we find useful 
to present here its full and consistent version. 



IV. CALCULATION SCHEME FOR M=3 



i.e. the transformation matrix for the phases is the in- 
verse transposed matrix for the populations. 



The triple-well system with interacting BEC and 
scheme of the cyclic STIRAP transport used in this study 
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are depicted in Fig. 2. As is seen from the plot a), the ad- 
joining wells are (weakly) coupled via the tunneling and 
the couplings Q12 and f2 2 3 play the role of the pump f2p 
and Stokes f^s Rabi frequencies in a familiar STIRAP 
scheme in Fig la). Further, the plot b) shows that, com- 
pared with other STIRAP applications for the transport 
of individual atoms [H| and BEC [H, [H HI , we imple- 
ment here three STIRAP transfers, |1) |3), |3) -> |2) 
and 1 2) — > |1), and every transfer follows the scheme a). 
In other words, these three transfers run via the inter- 
mediate states 1 2), |1), and |3) and use pairs of the cou- 
plings (012,^23), (^31,^12) and (S! 2 3,^3i), respectively. 
As was mentioned above, we need three STIRAP steps to 
produce the cyclic evolution of the system and generate 
a geometric phase. Besides, this allows to test fidelity of 
STIRAP transport in the chain of transfers. 

The similarity between the present scenario and typical 
STIRAP can be additionally justified by the comparison 
of their Hamiltonians. Following Fig. 2a), Eqs. (fl~5)) for 
M=3 can be written as 



ihip k (t) = £T(t)V*(t) 
with the Hamiltonian 



H(t) = 




Q P (t) 

E 2 + A 2 (t) a s (t) 
n s (t) E 3 + A 3 (t) 



(28) 



(29) 



where Ak{t) — ANh(t) is the non-linear interaction 
contribution and flp(t) = — fii 2 (i)/2 and f^s(i) = 
— ^23(0/2 are the pump and Stokes couplings. The non- 
linear terms Ak(t) are detrimental for adiabatic transfer 
within STIRAP [ljj. If to omit them, then dHJ) fully 
coincides with STIRAP Hamiltonian ©. 

In our study the time-dependent part of the coupling 
(TT4")) has the Gauss form 



f^.(t)=exp{-^L^} 



(30) 



where tkj and Tkj are centroid and width parameters. 
This form is smooth which is important for adiabaticity 
of the process. 

Using |14]) and (J30]) one may amend the STIRAP adi- 
abatic condition (fj]) . Following [36], |45j one gets 



Slip 



n 3 



« 1 or fi r » 



\d\ 



V2T 2 



(31) 



where Q T is the average amplitude of the pump and 
Stokes couplings during the overlap time r and d = 
ti2 — t23 is the relative pump-Stokes shift. Usually one 
may take il T w 0.5K and r w T. Then (j3Tj) is cast to 



Kt 2 > y/2\d\ 



(32) 



which means that STIRAP needs a strong coupling am- 
plitude K and/or a long overlap time r. It is easy to see 
that (|3"2"|) remains to be the same for the scaled dimen- 
sionless time. 



Since we use the mean-field approximation, the number 
of atoms in BEC should be sufficiently large to neglect the 
quantum corrections. In the present study we suppose 
that N > 10 4 [1, [nj. The total number of atoms is 
included to the parameter A in (fl6]l. 



V. RESULTS AND DISCUSSION 

1. Populations and phases 

Results of the calculations are depicted in Figs. 3-6. 
In all the figures time is dimensionless. The coupling 
parameters are T = 5.4 and d = -5. Then for K = 1 we 
get KT 2 « 30 and \2\d\ w 7 and so keep the adiabatic 
condition (j3"2l . 

In Fig. 3 the populations Ni(t) during three STIRAP 
steps are exhibited. The calculations are performed at 
initial conditions (t=0) iVi = 1, N2 = N3 = and 
0i = 4>2 = 03 = 0. The sequence of the pairs of Stokes 
(first) and pump (second) couplings is given in the panel 
a) while other plots demonstrate evolution of Ni(t) at 
different values of the ratio A and detuning A. 

As is seen from the panel b), all three STIRAP steps 
are robust and complete for A = A = 0, i.e. without in- 
teraction U and detuning. The minor peaks observed at 
t=20, 80 and 140, i.e. during the transfers, mean a weak 
temporary population of the intermediate wells, which 
hints that the transport is not fully adiabatic. However, 
what is important for our aims, the transport is com- 
plete. And this takes place even after 3 STIRAP steps 
and despite deviations (though minor) from adiabaticity. 

The next panels c)-e) show that switching on the inter- 
action worsens the transport. The damage is negligible 
for A = 0.1 and does not exceed 20% for A = ±0.2 but 
STIRAP ruins for larger interaction (not shown). In any 
case, STIRAP is robust under a modest interaction and 
this holds without any detuning. The comparison of the 
cases d) and e) shows that the result somewhat depends 
on the interaction sign. 

The remaining three panels f)-h) demonstrate role of 
detuning in the non- interacting condensate (A = 0). It is 
seen that a considerable detuning spoils the transport 
and the result slightly depends on the detuning sign. 
Our calculations generally confirm that a weak detun- 
ing is not harmful and, in accordance with [rJ I22I [23| , 
may be even useful to amend slightly adiabaticity of the 
process. However, the detuning is obviously not oblig- 
atory. Moreover, in real conditions, the adiabaticity is 
never fully kept. Nevertheless, the adiabatic transport 
schemes should work if diabatic perturbations are not 
strong, which is confirmed by our results. 

In Fig. 4 the dependence of the transport on the rela- 
tive shift d of the pump and Stokes couplings is demon- 
strated. The results are obtained for A = A = and 
r=4.36. Both counterintuitive (Stokes precedes pump, 
d < 0) and intuitive (pump precedes Stokes, d > 0) 
sequences of the couplings are covered. It is seen that 
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FIG. 3: STIRAP transport of BEC in the circular well chain, a): Sequence of the pairs of Stokes (first) and pump (second) 
couplings £lkj responsible for 1-2 (solid line), 2-3 (dash line) and 2-3 (dotted line) couplings. b)-h): Evolution of the populations 
Ni (solid line) N2 (dash line), and N3 (dotted line) at different values of the ratio A and detuning A. 
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FIG. 4: Population N3 at different relative shift d between 
pump and Stokes couplings. The maximal populations at 
counterintuitive (d < 0) and intuitive (d > 0) coupling order 
correspond to STIRAP and non-adiabatic transport. 



adiabatic. What is remarkable, there is no any transfer 
without the shift, i.e. for d = 0. Being adiabatic, STI- 
RAP transfer is less sensitive to the parameters of the 
process and so is more preferable than its non-adiabatic 
counterpart. This is partly confirmed by Fig. 4 where 
the left adiabatic peak is broader than the right non- 
adiabatic one, hence less sensitivity to the shift d. 

In Fig. 5 the phases fa and phase differences 6i are 
given for the cases with and without the interaction U . 
It is seen that the interaction and corresponding non- 
linear effects drastically change both fa and (9.;. This 
conclusion generally agrees with the observations for the 
oscillating BEC where the interaction also strongly 

affects the phases. So the interaction can in principle 
be implemented (via the Feshbach resonance) to control 
geometric phases generated during BEC transport. 



the best result (complete transport) takes place for the 
counterintuitive order, d = —3, pertinent for STIRAP. 
The intuitive order also leads to the appreciable popula- 
tion, N3 ~ 80 — 90% at d = 3, though this transfer is not 



A. Geometric phases 

Being coherent, BEC provides an interesting possibil- 
ity to generate and investigate various geometric (topo- 
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the oscillation and transport dynamics can be described 
within the same general formalism given in [46] ]. In the 
present study, we exploit its version (32j. There, for the 
cyclic evolution during the time interval t — [0,T], the 
geometric phase is determined as the difference between 
the total and the dynamical phases 



where 



lg=lt- Id 



7t = arg ty(0) • m) , 

Id = Im f dt$(t) ■ i(t)) 
Jo 



(33) 

(34) 
(35) 



FIG. 5: Phases fa (i=l,2,3, left panels) and phase differences 
9i (i=l,2, right panels) calculated without (A = 0, upper pan- 
els) and with (A = 0.2, bottom panels) interaction. In all the 
panels A = 0. For i=l,2,3 the solid, dash, and dotted lines are 
used. The phases in the upper and bottom panels correspond 
to the populations in Figs. 3b) and 3d), respectively. 
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FIG. 6: Total jt (solid line), dynamical 7d (dash line), and 
geometric 7 9 (dotted line) phases for the cases without (a) 
and with (b) interaction U. In both panels A = 0. In the 
panel a) the unconventional geometric phase 7 S — ajd with 
a 7^ — 1 is obtained. 



logical) phases 7 S [29|-[32|]. These phases are known to 
be mainly determined by the path topology (unlike the 
dynamical phases t„ which depend on the process rate 
as well), see e.g. [4a]. Because of this feature, 7 ff is less 
sensitive to parameters of the process and so can serve 
as a reliable carrier of information [? ]. During last 
decades implementation of geometric phases in the so 
called geometric quantum computation has become ac- 
tual [H, HH HH • In this connection using BEC transport 
for exploration of various j g could be indeed of a keen 
interest. 

The geometric phases in BEC dynamics have been al- 
ready studied elsewhere [2!|-[32j]. However, these studies 
concerned the oscillating BEC and, by our knowledge, 
still nothing was done for j g in BEC transport problems. 

In principle, geometrical and dynamical phases in 



(36) 



is the state vector consisting of the components © of the 
condensate. 

In Fig. 6 results of our calculations of 7t, 7<j, and 7 S 
for the three-step STIRAP transport are presented. Due 
to a cyclic adiabatic evolution, the Berry phase [47j is 
produced. The cases with and without interaction (cor- 
responding to the protocols of Figs. 3b), 3d) and 5) 
are considered. It is seen that in both cases we observe 
proportionality of the geometric and dynamical phases, 
7 g ~ ajd- However, for A = 0.2 there is a large mutual 
compensation of the phases (a « —1) thus giving 7t « 0. 
Instead for A = a w —2 (for t > 90) and we gain the 
so called unconventional geometric phase [34[ • This gives 
a chance to determine j g in interference experiments via 
measurement of jt- Indeed, if 7t ^ and we know a 
(which hopefully slightly depends on the parameters of 
the process), then we directly get j g . 

It worth noting that j g ^ only if the path has a 
non-trivial topology. Such topology becomes vivid if we 
reformulate the theory in terms of the sit(3) generators, 
in analogy to the su(2) quasi-spin operators treated e.g. 
in . In that two-mode case the BEC dynamics is 
reduced to the motion of the tip of the quasi-spin vector 
on the surface of the Bloch sphere. Hence the non-trivial 
topology of the path (e.g. as compared with motion on 
the plane). In such presentation, 7 9 is determined by the 
solid angle subtended by the closed curve on the sphere 
surface. The similar representation can be built for BEC 
transport in the triple-well trap as well. 

Altogether, BEC transport within STIRAP protocol 
seems to be a useful tool to generate different geometric 
phase. Such transport can be realized for a variety of the 
process parameters [22|, [23[ . So, a manifold of geometric 
phases can be produced. In this connection, it would be 
interesting to look for the STIRAP protocol leading to 
7t 7^ at 7^ = 0, i.e. for the conditions where only j g is 
produced. 
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VI. CONCLUSIONS 

The Stimulated Raman Adiabatic Passage (STIRAP) 
is applied to irreversible transport of the Bose-Einstein 
condensate (BEC) in the triple-well trap. The basic fea- 
tures of STIRAP are sketched and analogy between two- 
photon and tunneling STIRAP scenarios is discussed. 
The relevant formalism is presented and specified for the 
transport problem. 

The calculations are performed for the cyclic transport 
of BEC by using three successive STIRAP steps. It is 
shown that STIRAP indeed produces a robust and com- 
plete transport. Besides, it remains effective at modest 
interaction between BEC atoms and related non-linearity 
of the problem. As compared with the previous STIRAP 
studies [H, H3, HI] , w e demonstrate that detuning (trap 
asymmetry) is not obligatory and, at its large magni- 
tude, can be even detrimental (though small detuning 
can slightly amend adiabaticity of the process). 

Note that full adiabaticity of STIRAP can be hardly 
ensured in BEC transport since the transferred atoms 
must in any case pass the intermediate well thus disturb- 
ing the adiabatic following. In this connection, we do 
not pursue the perfect adiabaticity. Instead, we demon- 
strate that complete and robust transport can be realized 
even under its (though modest) distortion. Moreover, we 
show that effective transport can take place even at intu- 



itive sequence of partly overlapping couplings when the 
process is strictly non-adiabatic. Note that at zero in- 
teractions our results are relevant for the transport of 
individual atoms. 

For the first time, we demonstrate evolution of phases 
of BEC fractions in STIRAP transport and show that 
they strictly depend on the interaction. The correspond- 
ing dynamical and geometric phases are also computed. 
It is shown that at some interaction we gain the uncon- 
ventional topological phase which is proportional to its 
dynamical counterpart and both them produce a large to- 
tal phase. This finding may be used to determine uncon- 
ventional topological phases by measuring the total phase 
in interference experiments. Altogether, our study show 
that STIRAP transport can be a perspective tool for 
generation and exploration of various geometric phases 
which in turn are now of a keen interest for quantum 
computing [HHii]. 
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